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Abstract: 

The problem of solving multi-objective linear-programming problems, by assuming that 
the decision maker has imprecise goals for objective functions, is addressed. Several methods 
have been proposed in literature in order to obtain efficient solutions and then Pareto optimal 
solutions to multi-objective programming problems under imprecise environment. Recently 
Jimenez and Bilbao showed that a fuzzy efficient solution may not be Pareto optimal solution 
in the case that one of the fuzzy goals is fully achieved. Next Wu, Liu and Lur redefined the 
membership function of fuzzy set theory and proposed a new two phase method. In this paper 
we show that this redefined membership function has not been used during computation 
according to definition. An example is shows that if the redefined membership functions are 
used properly; it may fail to generate a solution. Therefore a new function that is strictly 
monotonic over entire set of real numbers as well as preserve an amazing characteristic of 
membership function is defined and used to find Pareto optimal solution of multi objective 
programming problem. Classical membership function is a special case of this newly defined 
function. Numerical examples illustrate our procedure. 

Keywords: Multi-objective programming; Decision analysis; Fuzzy mathematical 
programming; efficient solution; Pareto-optimal solution; Two-phase method; T- 
characteristic function. 



1. Introduction 

Let us consider the following multi-objective linear programming problem (MOLPP) 
with k objective functions z, (jc), i = 1. . . k as 

Min z(x) = (z 1 (x),z 2 (x),...,z i (x)) 

subject to xe X (1) 

HereX ={ie R" : Ax<b,x>0},b = (b l ,b 2 ,...,b m )<E R m } and A is a mxnmatrix. By assuming 
that the decision maker (DM) has imprecise aspiration levels for each objective function z t (x) 

, i = 1. . . k of model (1), several methods have been proposed in literature for characterizing 
Pareto optimal solutions to MOLPP (1). 

One such approach involves usage of fuzzy set theory [13, 14]. Here a MOLPP is 
converted into single objective optimization problem to attain fuzzy efficient solution [11, 
12] and then Pareto optimal solution by using the fuzziness of the DM's aspiration with 
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respect to the goals of imprecise objective functions (and constraints till they are symmetric!) 
[17]. Recently the concepts of evolutionary techniques such as particle swarm optimization 
and genetic algorithm are proposed to solve multi objective optimization problems. Suitable 
definition of continuity and differentiability of fuzzy-valued function have also been used to 
study multi-objective programming problems with fuzzy-valued objective functions. More 
solution concepts and techniques such as the KKT optimality conditions, the necessary and 
sufficient conditions of Pareto-optimality for imprecise multi-objective programming 
problems have also been discussed and tackled with. 

One such method to solve fuzzy MOLPP is two phase method [2,3]. Although the max- 
min operator is usually applied due to its easy computation, it does not guarantee to yield a 
fuzzy-efficient solution if more than one optimal solution exists [11]. The concept of two- 
phase procedure can also be employed to solve MOLPP with all fuzzy coefficients [10]. 
Sakawa contained rather comprehensive survey regarding interactive methods to MOLPP and 
fuzzy MOLPP [7, 8]. 

Recently Jimenez and Bilbao (2009) [6] showed that a fuzzy-efficient solution [11, 12] 
may not be a Pareto-optimal solution in the case that one of fuzzy goals is fully achieved. 
And they proposed a modified procedure which extended the two-phase approach of Guu and 
Wu (1997, 1999) [3, 4] and approach of Dubois and Fortemps (1999) [2] to attain a Pareto- 
optimal solution for model (1). 

But according to Yan-Kuen Wu et al (2015) [15] the proposed approach by Jimenez and 
Bilbao (2009) [6] cannot guarantee to be a general procedure to attain Pareto-optimal 
solution. And accordingly another two-phase approach is proposed to directly obtain a 
Pareto-optimal solution. According to them, the computation of this new two-phase approach 
is more simplified than Jimenez and Bilbao's modified procedure. In fact, in proposed 
approach by Wu et al, membership function of fuzzy set theory has been redefined [15]. 

The paper is organized as follows. In section 2, we identify that this redefined 
membership function has not been used according to its definition by Wu et al (2015) [15]. 
We also identify that a problem may become unnecessarily infeasible by keeping lower 
bound of membership function intact. In subsection 2.1., we discuss two examples and 
highlight that it is not always realistic to quantify or measure each type of imprecise 
information with in some bounded subset of real line. In subsection 2.2., a new function, call 
it T-characteristic function, and then T-efficient solution is introduced. In section 3, few 
lemmas are discussed related to Pareto optimal solution of an MOLPP under imprecise 
environment. Then in section 4, we propose an algorithm to find T-efficient and finally 
Pareto optimal solution to a MOLPP under imprecise environment. To further explain our 
approach, we solve numerical examples. The second example is same as suggested by 
Jimenez and Bilbao in 2009 [6]. The third one is same as suggested by Zimmermann in 1978 
[13]. Conclusions are drawn at last. 

2. Overview on membership function of fuzzy set theory 

In proposed approach by Wu et al in 2015 [15], membership function of fuzzy set theory 
is redefined. But their proposed mathematical models as well as numerical examples were not 
according to their definition of redefined membership function // ; (z ; (x)) . From mathematical 
model as well as from those two numerical examples, it is clear that only one part of 
redefined membership function // ; (z ; (x)) , where objective value does not exceed sum of goal 



and tolerance (for minimization type of objective functions J, was used. But these extra 
constraints were not added to set of existing constraints. Thus in problem formulation it 
happened that redefined membership function (z,. (*)) was strictly monotone over entire set 
of real numbers where as it was not so in definition. Later example 1 is given in this paper to 
show that those additional constraints may make a problem infeasible. 

Actually upper bound of membership function was removed by Wu et al (2015) but there 
still remained a lower bound at 0 [15]. In those examples by Wu et al, those additional 
constraints in definition of redefined membership function that objective value does not 
exceed sum of goal and tolerance (for minimization type of objective functions) were of little 
use [15]. But in some cases, if the lower bound of membership function remains intact, the 
problem may become infeasible. A solution is always better than no solution. "Half a loaf is 
better than none." A point to be noted. 

Thus the attention must be put on membership function, now. 

2.1. Membership function in fuzzy set 

Historically fuzzy set theory was developed from crisp set theory to quantify imprecise 
information. In crisp set, the belongingness of an element is represented by characteristic 
function^ : S — > {0,1} . An element belongs to a crisp subset A of universal set S if its value is 

1 and does not belong to the subset A if its value is zero. On the other hand, fuzzy subset A of 
universal set S can be characterized as a set of ordered pair of element* and its grade 

ju A (x) and is written as A = {< x,ju^(x) >: xe S] . Here the membership function //- is 

defined as//- : S — >[0,1], where S is universal set. According to Zadeh, the membership 

value ju A (x) indicates the degree of belongingness of the element x to the fuzzy set A. The 
membership values between 0 and 1 characterize members of S that belong to the fuzzy set 
A partially. 

Let us consider few examples. Consider fuzzy set A of all countries that have chance to 
win next Soccer World Cup in 2018. Most of us will agree that the chance of Russia is nil. So 

it can be safely stated that Russia is in A with membership value 0. Again chance of India's 

winning is also nil. So India is also in A with membership value 0. It means that Russia and 
India are at same level! But it is certain that even if Russia wins the 2018 Soccer World Cup, 
unfortunately India will certainly not be the winner. So in this case membership function of 
fuzzy set theory fails to explain! 

Again, suppose A is a fuzzy set representing good students. A student getting more than 
75% in examination may be considered a good student. So one scoring 76% has membership 
value 1 in fuzzy set A and another student scoring 98% will also have membership value 1 in 
fuzzy set A . Certainty, they are not same! 

Further analysis highlights that the cause lies in the definition of membership function 
itself. Membership function in fuzzy set theory is useful and provides satisficing result when 
extreme ends of imprecise information can be quantified within the boundary of zero and 
one. But it is not always correct or feasible to quantify or measure each type of imprecise 
information with in some bounded subset of real line. 



2.2. T-characteristic function, Pareto-optimal solution and T-efficient solution 

It is noted that an element of a fuzzy set can lie partially or never lie or must lie in that 
set. Idea of lying partially is well quantified or measured by membership function of fuzzy set 
theory. But it does not suit in case of must lying (membership value always one) and never 
lying (membership value always zero). 

Again in the case of decision making under fuzzy environment [1], one interesting and 
useful property of membership function is that higher value of membership function always 
gives better result of objective function. And there is not only membership function but also 
many such functions with this amazing characteristic. 

Therefore one new function, say T-characteristic function, that handles these underlying 
issues well, is defined. 

Definition 1. Let S denotes the universal set. Then T-characteristic function is denoted by 
T A (x) and is defined as 7^ : S — > R that assigns a real number T A (x) to each element xe S , 

where T A (x) represents the grade of membership of 5 in A. 

As for example, in decision making under imprecise environment [1], in case of 
maximization type of objective function zf, one may simply consider the T-characteristic 
function as z\ itself. And in case of minimization type of objective function Z2, one may 
consider T-characteristic function as (-) Z2. They satisfy the same property and can be used to 
convert a MOLPP into single objective optimization problem to obtain Pareto optimal 
solution. Next in imprecise environment, Pareto optimal solution and T-efficient solution are 
defined as follows: 

Definition 2. A decision plan x 0 e X is said to be a Pareto-optimal solution to the MOLP 
model (1) if there does not exist another ye X such that z ; (y)) < z ; (x 0 )for all i, i^j and 
Zj(y) < Zj(x 0 ) for at least one j. 

Definition 3. A decision plan x 0 e X is said to be a T-efficient solution to the MOLP model 
(1) if there does not exist anotherye X such that 7](z ( .(x 0 )) < r(z ; (y)) for all i, i^j and 
r. ( Z . (x 0 )) < Tjizjiy)) for at least one j. 

3. Obtaining T-efficient and Pareto-optimal solution 

Within the scope of the multi-objective decision making (MODM) theory, the Pareto- 
optimality of solution is a necessary condition in order to guarantee the rationality of a 
decision. Then using the method as discussed by Bellman and Zadeh, model (1) is equivalent 
to 

max (7; (z, (x)), T 2 (z 2 {x)),...,T k (z k (x))) 

si.xeX (la) 



And model (la) is equivalent to 



Min (7; ( Zl (x)), T 2 (z 2 (x)),...,T k (z k (x)))" 1 

s.t.xeX (lb) 
Using the min-max method by Bowman (1976) [16], this MOLPP (lb) may be rewritten 

as 

Min max (T t (z t (x)))" 1 , i = 1, 2.. .k 

i 

s.t.xeX (lc) 
Let v =max(T j (z i (x))y 1 J = l,2...k . Then model (lc) reduces to 

min v 
s.t. 

(T i (z i (xW l <v,i = l,2...k 

xe X, v unrestricted in sign (Id) 



If we put v = j/^,then (Id) reduces to 



max X 

s.t.T i (z i (x))>A,i = l,2,...,k 

x e X, X unrestricted in sign (2) 

The relationships between the optimal solution x* of the problem (2) and the Pareto 
optimal solution to the MOLPP (1) can be characterized by following lemmas. 

Lemma 1. If x*^X is unique optimal solution to the max-min problem (2), then x* is a 
Pareto optimal solution to the MOLPP (1). 

Proof. Let x* be unique optimal solution to problem (2). If possible let x* is not a Pareto 
optimal solution to problem (1). Then there exist at least one ye X such that z t (y)) ^ z,(x*) 
for all i, i j and Zj(y)< z y (x*)for at least one j. Since each r(z,)is strictly monotonic 
decreasing function for minimizing type of objective functions z ; , i= 1, 2. . .k, we get 

7](z,.(x*)) <7;(z,.(y)) for all i =1, 2,...,k , i * j 
and r.(z 7 (x*))<r.(z 7 (y)). 

It contradicts that x* is unique optimal solution to problem (2). Therefore our assumption 
is wrong. Hence x* is Pareto optimal solution to MOLPP ■ 

Lemma 2. If x* E X is a Pareto optimal solution of the MOLPP, then x* is an optimal 
solution of the problem (2) for some fi, i=l, 2...k. 

Proof. Let x* is Pareto optimal solution to MOLPP. If possible let x* is not an optimum 
solution to the minimax problem for any 7^., i = 1, 2...k. 



Then for T , i-l,...k, there exists xe X such that T j (z j (x)) > T j (Zj (x*)) = A(x*), for some 
j and T t { Zi (x)) >T t { Zi (x*)) = X * (x*), i = \,...k,i* j. 

But 7^ is strictly decreasing function for each i = l...k. Hence it gives 
Zj(x) < z f {x*) for some j and z ( .(x) < z t {x*),i - \,2,...k,i ^ It contradicts that x* is Pareto 
optimal solution to the MOLPP. Therefore our assumption is wrong. Hence x* is an optimal 
solution to the minimax problem (2) for?] , i =1, 2. .. km 

From lemma 1 and lemma 2, if the uniqueness of the optimal solution x* for the minimax 
problem (2) is not guaranteed, it is necessary to perform the Pareto optimality test of x*. 

The Pareto optimality test for x* can be performed by solving the following linear 
programming problem with the decision variables x = (x 1 ,x 2 ,...x n f and e = (e l ,e 2 ,...e k ) T 
[2,3]- 

k 

max^e ; 

;=1 

T i (z i (x))-e i >T i (z i (x*)),i = l,2,...,k 

xe X,£ i >0,i = \,2,...,k. (3) 

Lemma 3. Let x and £ be optimal solution to problem (3). Then 

I. If all £ i - 0 , then x* is Pareto optimal solution of the MOLPP. 
II. If at least one £ i > 0 , then x* is not Pareto optimal solution of the MOLPP. Instead 
of x*, x is Pareto optimal solution to the MOLPP. 

Proof of this lemma is already recorded in literature [3, 4]. It may be noted that several 
authors have proposed other procedures to test Pareto optimality and can be used in place of 
model (3). 

4. A general framework for solving multi-objective programming problems under 
imprecise environment 

The above ideas can be integrated into a general framework or algorithm to obtain a 
Pareto-optimal solution to the MOLPP which has the additional property of being T-efficient 
for chosen T-characteristic functions. The steps of the proposed algorithm can be synthesized 
as follows [9, 10] 

Step 1: Define suitable T-characteristic function for each imprecise objective function in 
such a way that higher value of T-characteristic function always gives better value of 
objective function. Well-defined T-characteristic functions preserve this characteristic and 
gives finite value always. 

Step 2: Construct the max-min problem as model (2). 

Step 3: Solve the max-min problem and obtain the T-efficient solution x* of it. 

Step 4: To test whether this T-efficient solution x* is Pareto optimal, solve model (3). 

Step 5: Let x and £ be optimal solution of model (3) in step 4. Then two cases may arise: 



a) If e i - 0 V i , then x* is Pareto optimal solution of the MOLPP. 

b) If e > 0 for at least one j , then x* is not Pareto optimal solution of the 
MOLPP. Instead of x*, x is Pareto optimal solution to the MOLPP. 

Step 6: This solution is T-efficient as well as Pareto-optimal. The algorithm is complete. 

4.1. Numerical examples 

To check the usefulness of proposed approach, consider these numerical examples. The 
first example is chosen to show that the addition of constraints, necessary to make redefined 
membership function consistent with its definition, as proposed by Wu et al [15], may make 
problem infeasible. The second one is the same as suggested by Jimenez and Bilbao in 2009 
[6]. The third one is the same as suggested by Zimmermann in 1978 [13]. 

4.1.1. Example 1 

Let us consider a fuzzy multi objective linear programming problem as 

fuzzy max 5x l + 5x 2 
fuzzy min 5x x + x 2 
fuzzy max 3xj-8.2x 2 
s.t. 

5x l +lx 2 < 12, 9x l + x 2 < 10, - 5x 1 + x 2 < 3, 
x l , x 2 > 0 

4.1.1.1. Using redefined membership function 

Suppose that the DM specifies goal and tolerance for each of the imprecise objective 
functions as follows: 



Objective function 


Goal 


Goal plus Tolerance 


Fuzzy Maximum 


10 


9 


Fuzzy Minimum 


0 


1 


Fuzzy Maximum 


3 


-7 



Table 1. Goal and tolerance of fuzzy objective functions 

It must be noted that these values of goal and goal plus tolerance of each of the objective 
functions are within their individual maximum and minimum values. Based on this 
information, as Wu et al [15], redefined membership functions are defined as 



ju 1 (z 1 (x)) = 



5x 1 +5x 2 9 9 
1 1 2 

0 ,5x,+5x 2 <9 



// 2 (z 2 (x)) = - 



l-(5jtj + x 2 ) 

_ 



,5jCj + x 2 < 1 
, 5Xj + x 2 > 1 



// 3 (z 3 (x)) = 



3Xl S2Xl 3 ,3x t -8.2x 2 >-7 



10 
0 



,3xj -8.2x 2 < -7 



The max-min problem, using the redefined membership function properly (according to 
their definition), is as follows: 

max A 
s.t. 

5x +5x, —9^ r> l-(5x,+x 9 )^ , 3x,-8.2x 7 -3^ - 

— > A, — > A, — > A, 

1 1 10 

5x l +5x 2 >9, 5xj + x 2 <1, 3x r 8.2x 2 >-7, 

5xj+7x 2 <12, 9x!+x 2 <10, -5x,+3x 2 <3, 
x 1 ,x 2 > 0, A is unrestricted in sign 

Using Lingo 15.0.32 to solve this problem, we find that the problem has no feasible 
solution. 

4.1.1.2. Using T-characteristic function 

Now the T-characteristic functions for each of the objective are formed, keeping in mind 
that the higher value of T-characteristic function gives better value of corresponding 
objective function in each case as 

r 1 (z 1 (x)) = 5x ' + ^ 2 " 9 ,x 1 ,x 2 >o, 

r 2 (z 2 (x)) = 1 ~ (5 ^ +X2) ,x 1 ,x 2 >0, 
r 3 (z 3 (x)) = ^|^,x 1 ,x 2 >0 



If these T-characteristic functions are used instead of membership functions, the problem 
in model (2) may be formulated as 

max A 
s.t. 

5x,+5x 9 -9^ , l-(5x, +x 9 *> , 3x -8.2x 9 -3 ^ , 

— > A, — > A, — > A, 

1 1 10 

5x!+7x 2 <12, 9x 1 +x 2 <10, -5x^3x2 < 3, 
x p x 2 > 0, A is unrestricted in sign 



Using Lingo 15.0.32 to solve this problem, the optimal solution is obtained as 

X* = -1.333333, x* = 0.2, x* = 1.333333, z, * (x*) = 7.665, z 2 * (x*) = 2.333, z 3 * (x*) = -10.33 1, 
T x * ( Zl * (x*)) = -1 .333333, T 2 * (z 2 * (x*)) = -1.333333, J 3 * (z 3 * (**)) = -1.333333. 

It is observed that values of T-characteristic functions are negative. To test the Pareto 
optimality of this solution, construct model (3) for this problem as follows 



_ 3 

max 

!=1 



e, > -1.333333, 
e 2 > -1.333333, 
-e 3 > -1.333333, 



subject to 

- 

l-(5x, +x 2 ) 
_ — 

3xj -8.2x 2 -3 
10 

5xj+7x 2 <12, 9x!+x 2 <10, -5x!+3x 2 <3, 
x,,f,.>0, i = l,2,3. 

Let xandfbe optimal solution to this problem. Using Lingo 15.0.32 to solve this 
problem, the optimal solution is obtained as follows: 

e x =e 2 =e 3 =0,x\ =0.2,x 2 =1.333333, z l (x) = 7.665, 1 2 (x) = 2.333, z 3 (x) = -10.331. 



Here all e t - 0. Hence the solution x* obtained in previous stage is itself Pareto optimal 

solution. It is observed that the proper usage of definition of redefined membership function, 
as proposed by Wu et al (2015) [15], makes this problem infeasible, whereas Pareto optimal 
solution is found by using proposed T-characteristic function. 

4.1.2. Example 2 

To illustrate the proposed method further, consider the numerical example suggested by 
Jimenez and Bilbao (2009) [6] as follows: 

Minimize z ; (x) = 3x t + 3x 2 + 3x 3 
Minimize z 2 (x) = 2xj+ x 2 + 2x 3 
Minimize z 3 (x) = 4x ; + 4x 2 + 2x 3 
subject to 4Xj + 2x 2 + 4x 3 > 18 
x t >l, x 2 >0, 0<x 3 <3 

4.1.2.1. T-characteristic function similar to membership function of fuzzy set theory 

The proposed approach to solve this MOLPP with fuzzy goals g x = 21, g 2 =S,g 1 =13 and 
tolerances t x -3,t 2 =2,t 3 =2 is employed. Based on proposed definition, the T-characteristic 
functions (partly similar to membership function) may be defined as follows: 



, 24-(3x l +3x,+3x 3 ) ^ . XN 10-(2x+x,+2x,) 

r 3 ( z, (.»-)) = 1 5 ~ (4 *' + * x i +2Xl) , Xl ,x,,x,>o 

Then the model (2) can be formulated as below: 
max A 

subject to the constraints 

24 -(3x 1 +3x 2 +3x 3 ) 10-(2x 1 + x 2 + 2x 3 ) 15-(4x 1 +4x 2 +2x 3 ) 

> A, > A, > A, 

111 

4xl+2x 2 +4x 3 >18, X;>1, 0<x 3 <3, 
x 2 > 0, A is unrestricted in sign 

Using Lingo 15.0.32 to solve this problem, the optimal solution is obtained as 

A* = 0.5, x* = 2.5, x 2 * = 0, x 3 * = 2, z x * (x*) = 13.5, z 2 * (x*) = 9, z 3 * (x*) = 14, 
T x * (z, * (x*)) = 3.5, T 2 * (z 2 * (x*)) = 0.5, r 3 * (z 3 * (x*)) = 0.5. 

To test the Pareto optimality of this solution, construct model (3) of this problem as 
follows 

3 

max V e i 

i=\ 

subject to the constraints 



24-(3x 1 +3x 2 +3x 3 ) 
3 

10-(2x 1 +x 2 + 2x 3 ) 



2 



-e l >3.5, 
-e >0.5 



15-(4x 1 +4x 2 +2x 3 ) £ ^ Q5 



2 

4x 1 +2x 2 +4x 3 >18, Xj>1, 0<x 3 <3, 
x 2 ,£-,. > 0,i =1,2,3 

Let xandfbe optimal solution to this problem. Using Lingo 15.0.32 to solve this 
problem, we obtain the optimal solution as 

e x =e 2 =0,e 3 =1,X[ =1.5,x 2 = 0,x 3 = 3,z 1 (x 1 ) = 13.5,z 2 (x 2 ) = 9,z 3 (x 3 ) = 12. 

It is observed thatz 3 -< z 3 *■ And according to lemma 3, x is the Pareto optimal solution. It 

is further noted that this Pareto optimal solution x is same as solution obtained by Wu et al 
[15]. 

4.1.2.1.1. Another T-characteristic function 



Again the T-characteristic functions may be defined for minimization type of objective 
function as 



Tj (Zj (x)) = -(3x 1 + 3x 2 + 3x 3 ), Xj , x 2 , x 3 > 0 



T 2 (z 2 (x)) = -(2xj + x 2 + 2x 3 ), Xj , x 2 , x 3 > 0 
T 3 (z 3 (x)) = -(4xj + 4x 2 + 2x 3 ), Xj , x 2 , x 3 > 0 
Based on these, model (2) can be formulated as 
max X 

subject to the constraints 

-(3xj + 3x 2 + 3x 3 ) > X, - (2xj + x 2 + 2x 3 ) > X, - (4x l + 4x 2 + 2x 3 ) > X 
4x!+2x 2 +4x 3 >18, x t >l, 0<x 3 <3, 
x 2 >0,X is unrestricted in sign 

Using Lingo 15.0.32 to solve this problem, the optimal solution is obtained as 

X* = -13.5,V = 1.5,x 2 * = 0, x 3 * = 3, z* = 13.5, z* = 9, z 3 * = 12, 
T x * (z, * (x*)) = -13.5, T 2 * (z 2 * (x*)) = -9, T 3 * (z 3 * (x*)) = -12. 

To test the Pareto optimality of this problem, construct model (3) of this problem as 
follows 

3 

max^V, 

;=1 

subject to the constraints 
-(3x; +3x 2 +3x 3 )-£"j > —13.5 
-(2x l + x 2 + 2x 3 ) - e 2 > -9 
-(4xj +4x 2 +2x 3 )-£- 3 >-12 
4xj+2x 2 +4x 3 > 18, x t >l, 0<x 3 <3, 
x 2 ,e i >0,i = 1,2,3. 

Let xandfbe optimal solution to this problem. Using Lingo 15.0.32 to solve this 
problem, obtain the optimal solution as 

£^ — £ 2 — ^3 — — 1*5} ^2 — 0} jc^ — 3 5 — 1 3 »5j ^2 — 9 5 — 1 2) 
TOCO) = -13.5,f 2 (z 2 (x)) = -9,f 3 (z 3 (x)) = -12. 

Here all e i - 0. Hence the solution obtained in previous stage is itself Pareto optimal 
solution. 

4.1.3. Example 3 

A company manufactures two products 1 and 2 on given capacities. Product 1 yields of 
profit of 2 $ per piece and product 2 of 1 $ per piece. While product 2 can be exported, 
yielding a revenue of 2 $ per piece in foreign countries, product 1 needs imported raw 
materials of 1 $ per piece. Two goals are established: (a) Profit maximization and (b) 



Maximum improvement of the balance of trade, i.e. maximum difference of exports minus 
imports. 

Then the problem can be modelled as 

max imize z l (x) = -x t + 2x 2 
maximize z 2 (x) = 2x t +x 2 ; 
subject to 

-x 1 +3x 2 <21, Xj+3x 2 <27, 4xj+3x 2 <45 
3xj +x 2 < 30, Xj,x 2 > 0 

Suppose we construct T-characteristic function T(z(x)) of objective function z(x) of 
maximization type (partly similar to membership function of fuzzy set theory, as used by 
Zimmermann) as 

z +3 z —1 

T l (z l (x)) = -^-,x l ,x 2 >0, T 2 (z 2 (x)) = ^-,x l ,x 2 >0. 

Based on these functions, model (2) can be formulated as follows: 

max A 
subject to 

-x t + 2x 2 +3 ^ 2x, + x 2 -7 ^ 
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-x 1 +3x 2 <21, x t +3x 2 <27, 

4x 1 +3x 2 <45, 3x^X2^30, 
x p x 2 > 0,A is unrestricted in sign 

Using Lingo 15.0.32 to solve this problem, the optimal solution is obtained as 

X* = 0J4,x* = 5.03, x* = 7.32, z,* = 9.61, z 2 * = 17.39, 
T x * (z, * (x*)) = 0.74, T 2 * (z 2 * (**)) = 0.74. 

To test the Pareto optimality of this problem, construct model (3) of this problem as 
follows 

2 

max ^ e i 

i=\ 

subject to 

-x t + 2x 2 +3 > ^ 2x 1 + x 2 -7 > 

17 ' 14 

-Xj+3x 2 < 21, x t +3x 2 < 27, 

4x t +3x 2 <45, 3x t +x 2 <30, 
x„£ >0,* = 1,2. 



Let x and e be optimal solution to this problem. Using Lingo 15.0.32 to solve this 
problem, obtain the optimal solution as 

e 1= e 2 =0,x, = 5.03, x 2 = 7.32,^ =9.61, z 2 =17.39, 
?i(z 1 (x)) = 0.74,f 2 (z 2 (x)) = 0.74. 

Here all £ i - 0. Hence the solution obtained in previous stage is itself Pareto optimal 

solution. The maximum degree of "overall satisfaction" (0.74) is achieved for the solution x 
= (5.03, 7.32) T . And interestingly this result is exactly matching with result obtained by 
Zimmermann in 1978 [13]. 

5. Conclusions 

In this paper, we identify that the redefined membership function, proposed by Wu et al 
(2015) [15], was not used according to definition in their paper. It is not strictly monotonic 
over its entire domain in the definition (lower bound zero) but is so in problem formulation as 
well as in those two numerical examples (constraint showing lower bound is not present). 

In this paper we further identify that although this redefined membership function 
correctly identifies the best solution among the good; it fails to identify bad solution among 
the worst! As shown in ex. 1 of our paper, proper usage of redefined membership function 
may make the MOLPP infeasible. 

Moreover in real life situations, DM may not be able to identify the goal or tolerance or 
both; and in some cases, a DM may not exist at all. In these scenarios, it is problematic or 
even impossible to construct fuzzy membership function or that redefined membership 
function. 

Further we note that the primary objective is to identify a function such that higher value 
of function will generate better value for objective. Hence in decision making under 
uncertainty, we propose to define and use T-characteristic function covering the entire set of 
possible solutions in place of membership function. And this is a general procedure to obtain 
a Pareto-optimal solution that has the additional property of being T-efficient. 

And in fact the usage of a well-defined T-characteristic function also removes the need to 
collect information from DM and saves precious time of both. 

From this we also conclude that, under these circumstances, the issue of getting no 
solution even by using fuzzy set theory yielded from several published methods seems 
worthwhile or even necessary to reconsider. The same is also true for nonlinear problems of 
this type. 
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